Effects of magnetism and electric field on the energy gap of bilayer graphene 

nanoflakes 
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We study the effect of magnetism and perpendicular external electric field strengths on the energy 
gap of length confined bilayer graphene nanoribbons (or nanoflakes) as a function of ribbon width 
and length using a first principles density functional electronic structure method and a semi-local 
exchange-correlation approximation. We assume AB (Bernal) bilayer stacking and consider both 
armchair and zigzag edges, and for each edge type, we consider the two edge alignments, namely, 
a and /3 edge alignment. For the armchair nanoflakes we identify three distinct classes of bilayer 
energy gaps, determined by the number of carbon chains in the width direction (N — 3p, 3p+l and 
3p+2, p is an integer), and the gaps decrease with increasing width except for class 3p+2 armchair 
nanoribbons. Metallic-like behavior seen in armchair bilayer nanoribbons are found to be absent in 
armchair nanoflakes. Class 3p+2 armchair nanoflakes show significant length dependence. We find 
that the gaps decrease with the applied electric fields due to large intrinsic gap of the nanoflake. 
The existence of a critical gap with respect to the applied field, therefore, is not predicted by our 
calculations. Magnetism between the layers plays a major role in enhancing the gap values resulting 
from the geometrical confinement, hinting at an interplay of magnetism and geometrical confinement 
in finite size bilayer graphene. 

PACS numbers: 71.15.Mb, 81.05.Uw, 75.75.+a 
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I. INTRODUCTION 

Bilayer graphene continues to attract attention from 
the condensed matter physics community as well as from 
industry because of the prediction and subsequent exper- 
imental realization of electric field opening of the band 
gapji the theoretical predictions of a novel electronic 
phase and the theoretical realization of a novel charge 
based electronic switch^ The modulation of bilayer gaps 
by a geometrical confinement (in particular by widths), 
interlayer magnetism and edge alignments 3 ' offers excit- 
ing opportunities for designing novel electronic and opti- 
cal devices. 

Recent advancement in the growth of large area 
graphene films on metal substrates, 5 deposition of a di- 
electric layer on the graphene film^ and making contacts 
to graphene for transport studies 7 motivates the theo- 
retical understanding of the intrinsic (edges and their 
alignments, edge disorders, layer stackings) as well as 
extrinsic (substrates, contacts, dielectrics, adatoms, de- 
fects) perturbations to the ideal properties of graphene 
and graphene sheets. In transport measurements, the 
edges of the graphene flakes which are in electrical con- 
tact with a piece of metal or an another piece of graphene 
play a decisive role in stabilizing an interface potential 
barrier for flow of electrons across it and this potential 
barrier depends on the intrinsic gap of the graphene flake. 

Theoretical studies of isolated finite size graphene 
flakes are, therefore, critical to understand the tunability 
of the gap value resulting from the edge confinements. 
Moreover, magnetism which is absent in bulk graphene 
bilayer can appear because of edges in the finite size 
graphene flakes. In addition, in bilayer systems, one can 



apply an external electric field to tune the energy band 
gap»i It is, therefore, desirable to understand the effects 
of magnetism and external electric fields on the energy 
gaps of geometrically confined isolated graphene sheets. 




armchair 



FIG. 1: (Color online) Schematic illustration of the two edge 
alignments we consider in bilayer graphene nanoflakes. (a) 
a- alignment and (b) /3- alignment. 

In this article, we use a density functional the- 
ory based first principles method 8 and a version of 
semi-local exchange correlation potential called gradi- 
ent approximation^ to explore the dependence of the 
energy gaps of bilayer graphene on geometrical confine- 
ments (length as well as width), interlayer magnetism 
and the external electric fields. We assume Bernal 
(or AB) stacked bilayers with both edge types (arm- 
chair and zigzag) and for each edge type, we consider 
two edge alignments (Fig. [1]), denoted as a- and /3- 
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alignments. Such length confined nanoribbons will be 
called "nano flakes" throughout this article. In the follow- 
ing discussions, we will refer to gradient approximation 
as GGA. 

We are not aware of any density functional theoret- 
ical calculations focussing on interplay of energy gap, 
magnetism and external electric field in bilayer graphene 
nanoflakcs in the literature. Only recently, the energy 
gap of monolayer nanoflakes was studied using a density 
functional based approach^ 

Here, we predict interplay of magnetism and geometri- 
cal confinement on the energy gap, the existence of three 
distinct classes of gaps in armchair bilayer nanoflakcs, 
and external electric field decreasing the energy gaps, in 
addition to capturing some of the predictions made for 
monolayer graphene nanoflakes. We note that our pre- 
dicted gap values are useful in optical experiments in 
which the possibility of exciton formation can be ruled 
out. The presence of bound excitons in an absorption 
spectra, however, can significantly change the intrinsic 
gaps of the graphene flakes. 

Our paper is organized as follows. We first summarize 
the density functional theory (DFT) electronic structure 
calculations that we have performed, commenting on the 
motivation for choosing GGA and the different combina- 
tions of lengths and widths in section II. In section III, 
we present the results that we have obtained for bilayer 
nanoflake gaps, focusing most extensively on the inter- 
play between width and length, edge magnetism, and 
the external electric field between the layers. Finally we 
summarize our results and present our conclusions. 



II. DENSITY FUNCTIONAL THEORY 
CALCULATIONS 

Our electronic structure calculations are per- 
formed with plane wave basis sets and ultrasoft 
pscudopotentialsAi implemented in a density functional 
theory based electronic structure method^ We use 
GGA- as our previous study 3 suggested the inadequacy 
of the local density approximation (LDA)i2 in capturing 
the magnetic nature of the ribbons. We note that, in 
the case of a-aligned zigzag nanoribbons, the magnetic 
ground state was found to be sensitive to the nature of 
the exchange-correlation potential^ Also the choice of 
GGA enabled us to compare our bilayer nanoflake results 
with those of the monolayer graphene nanoflakcs. 10 

We fixed the interlayer separation to 0.335 nm between 
the nanoflake sheets as GGA is found to highly overesti- 
mate the interlayer separation of the graphene sheets^ 
This is due to the absence of non-local or van der Waals' 
interactions in the LDA, GGA, Perdew-Burke-Ernzerhof 
(PBE) 15 and the hybrid versions of exchange-correction 
potential in DFT. We placed the nanoflakes in a super- 
cell set-up by introducing vacuum regions of 1 nm along 
all three superecell directions to avoid the intercell in- 
teractions. We used 2x2x1 k-point mesh in the full 



Brillouin Zone (BZ) and 476 eV kinetic energy cut-offs. 
The convergence of total energy was carefully tested with 
larger energy cut-off, vacuum size and the k-point mesh. 
Since we are considering finite size systems in a super- 
cell set-up, the energy gaps are computed at the T point 
of the irreducible BZ. To study the behavior of the gaps 
in the presence of external electric field, we used several 
values of electric field, perpendicular to the layers, up to 
a maximum value which is close to the Si02 dielectric 
breakdown field of 1 V/nm^ 

The in-plane er-orbitals at the edges were saturated 
with hydrogen atoms (with the C-H distance chosen to 
be the C-H bondlength in the CH4 molecule). Since 
the magnetism arises purely due to the edge geometry 
and localized carbon 7r-orbital, choosing a different C-H 
bondlength will not alter the interlayer magnetic order. 
However, we note that different edge functionalizations, 
other than hydrogens, can alter the electronic structure 
of the bilayer nanoflakes. Our predictions are therefore 
most relevant for nanoflakes cut in a hydrogen environ- 
ment without any significant edge disorder^ Moreover, 
we do not consider edge disorder and the role of different 
substrates in our calculations. 

In order to find the nature of interlayer magnetic order 
in bilayer nanoflakcs, we performed three distinct sets of 
calculations: a nonmagnetic, a parallel and an antipar- 
allel arrangement of magnetic moments on the carbon 
atoms between the layers. The intra-edge magnetic or- 
der within a layer was assumed to be antiparallel as var- 
ious theoretical reports suggests^ To find the magnetic 
ground state for bilayer armchair and zigzag nanoflakes, 
we chose a few representative systems. We used both 
narrow and large width nanoflakes with two edge types 
and for each edge type, two edge alignments (Fig. [1]). 

We find that in all cases, the nanoflake with interlayer 
antiparallel (or antiferromagnetic) magnetic moment ar- 
rangement is energetically favorable*^ This ground state 
is lower, by about 0.4 eV, than the interlayer parallel ar- 
rangement of the moments, and the nonmagnetic energy 
barrier is about 1 eV from the ground state of the in- 
terlayer antiferromagnetic order. Therefore, all the gap 
calculations in this article are done with antiferromag- 
netic order between the layers*^ 

We note that for both edge types, nanoflakes with 
a-alignment is found to be energetically favorable over 
those with f3- alignments and for the same width and 
length, /3-aligned edges affect the energy gaps more 
strongly compared to their non-magnetic counterparts 
than the a-aligned edges. For sake of completeness, we 
will consider both edge alignments for each edge type in 
this paper. 



III. GRAPHENE BILAYER NANOFLAKE GAPS 

We now present our results for the width, length, mag- 
netism, and external electric field dependence of the bi- 
layer gaps in nanoflakes with both edge types and two 
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edge alignments. First we discuss the width depen- 
dence of both armchair and zigzag nanoflakes with a fixed 
lengths of L w 3.2 nm and 3.3 nm respectively. Then we 
vary the lengths up to 6.1 nm and discuss the dependence 
of the bilayer gaps on the length. It should be noted 
that an armchair nanoflake with the width/length ratio 
( W/L) can equivalently be defined as a zigzag nanoflake 
with the aspect ratio (L/W). We will close this section 
with the external electric field effects on the gaps with 
fixed as well as variable lengths. 



Armchair and zigzag bilayers with fixed length 
confinement 



In this section, we discuss the width dependence of 
gaps in bilayer graphene nanoflakes with both edge align- 
ments for a fixed length. Figure [2] shows variations of the 
gap versus the nanoflake width for a-aligned armchair 
nanoflakes. For comparison, we also show the gaps for 
bilayer armchair nanoribbons. We predict three classes 
of armchair gaps which we label, according to the widely 
accepted notations in the literature, by N = 3p, 3p+l 
and 3p+2, where p is an integer. As expected on the 
basis of previous workj* the width dependence is quite 
smooth within the three nanoflake classes, which are dis- 
tinguished by the number N of carbon chains across the 
nanoflake mod 3. However, all classes show a semicon- 
ducting behavior. This is in contrast with the class 3p+2 
in armchair nanoribbons which exhibited a tendency to- 
wards metallic behavior whereas the class 3p+2 in arm- 
chair nanoflakes remain semiconducting due to the mag- 
netic zigzag units present along the width direction. 

For the class 3p, magnetism between the bilayers sig- 
nificantly enhances the gaps compared not only to non- 
magnetic bilayer nanoflakes but also to corresponding 
nanoribbons, as seen in Fig. [U[a). But this is not im- 
mediately apparent for the class 3p+l until the width 
reaches about 3 nm, where the nanoflake gap tends to 
exceed the nanoribbon gap. To check whether for larger 
lengths and widths, magnetic nanoflake gaps exceed cor- 
responding nanoribbon gaps for 3p+l, we considered a 
longer ribbon (L « 4.2 nm) and a width near 4 nm and 
found that indeed the bilayer nanoflake gap is larger than 
the nanoribbon gap (0.27 eV versus 0.12 eV). This indi- 
cates an interplay of interlayer magnetism and geometri- 
cal confinement in bilayer armchair nanoflakes because 
the number of magnetic zigzag units increases as the 
width increases, which slows down the decrease of the 
armchair nanoflake energy gaps compared to nonmag- 
netic armchair nanoribbon energy gaps. The behavior of 
the class 3p+2 nanoflakes can be understood similarly. 
For this class, the length confinement resulted in larger 
gaps compared to the nanoribbons and the magnetism 
seems to be contributing strongly in enhancing the gaps 
for wider nanoflakes, which is opposite to the behavior of 
wider nanoflakes in two other classes 3p and 3p+l. Such 
increase of class 3p+2 gaps with the width was also seen 
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FIG. 2: (Color online) Variation of the energy gap with width 
of a-aligned bilayer armchair nanoflakes for a fixed length L 
~ 3.2 nm (open circles). Three classes of the nanoflakes are 
clearly seen in (a) 3p, (b) 3p+l and (c) 3p+2 where p is 
an integer specifying the number N of carbon chains along 
the width direction. Here p =1,2,- ■ ■ ,8, which translate to 
nanoflakes with widths less than but close to 3 nm. The 
gaps of the magnetic nanoflakes (open circles) are compared 
with those of the nonmagnetic nanoflakes (solid circles) and 
three classes of nonmagnetic nanoribbons (open triangles). It 
should be noted that armchair nanoribbons are all nonmag- 
netic and all the gaps are obtained using the same semi-local 
GGA potential. 



in monolayer nanoflakes 

DFT again predicts three classes of gaps for /3-aligned 
nanoflakes (Fig. Compared to the gaps in a-aligned 
nanoflakes, the gaps are consistently larger for the same 
width. Note that in bilayer nanoribbons^ the two types 
of edge alignments provide distinct electronic structures: 
magnetism in /3-aligned ribbons were more significant 
than in cv-aligned ribbons due to a dispersionless non- 
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FIG. 3: (Color online) Same as Fig. 2 but for /3-aligned 
bilayer armchair nanofiakes. 



magnetic band at the Fermi level spanning ^ of the dis- 
tance from the edge of the one dimensional BZ although 
inclusion of magnetic order opened up a gap in the energy 
spectra of both types. The two nanoflake classes 3p and 
3p+l show similar behavior as in a-aligned nanofiakes, 
again with a hint of order change in class 3p+l. In 
class 3p+2 nanofiakes, the gaps increase with the increase 
of widths for magnetic nanofiakes while nonmagnetic 
nanofiakes show behavior similar to the nanoribbons. For 
wider magnetic nanofiakes of W>2A nm, however, the 
energy gaps begin to decrease as width increases, again 
indicating the complex interplay of magnetism and geo- 
metrical confinement. 

Figure HI shows variations of the gap with the widths in 
a- and /3-aligned zigzag nanofiakes. For comparison, gaps 
of magnetic nanoribbons and nonmagnetic nanofiakes are 
also shown. Interlayer magnetism seems to enhance the 
gap values compared to the gaps resulting from only ge- 
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FIG. 4: (Color online) Variations of the energy gap of the 
(a) a-aligned and (b) /3-aligned bilayer zigzag nanofiakes with 
a fixed length L ft. 3.3 nm. For comparison, the gap val- 
ues of nonmagnetic nanofiakes and magnetic nanoribbons are 
also shown. All the gap calculations are done with the same 
semi-local GGA potential and for magnetic calculations an 
interlayer antiferromagnetic order was considered. 



ometrical confinement (i.e. nonmagnetic gaps), but in 
contrast to the armchair nanofiakes, the energy gaps of 
zigzag nanofiakes follow the same trend as corresponding 
zigzag nanoribbons. This may be due to the fact that 
with increasing width in armchair nanofiakes, the num- 
ber of (magnetic) zigzag units increases whereas in zigzag 
nanofiakes, increasing the width increases the number of 
(nonmagnetic) armchair units. 



B. Length dependence of the gaps 

We address the length dependence of the nanoflake 
gaps in this section. The metallic and non-metallic 
nanoribbons are expected to respond differently to 
changes in length as the authors of Ref. 10 show. There- 
fore, we considered length confined metallic (class 3p+2) 
and non- metallic (class 3p+l) bilayer armchair nanorib- 
bons as well as the non-metallic bilayer zigzag nanorib- 
bons for the study of length dependence of the gaps. We 
denote these, respectively, as metallic and non-metallic 
nanofiakes in Fig. [5] Both edge alignments were used and 
both magnetic and non-mangetic nanofiakes are consid- 
ered. We considered representative widths corresponding 
to N=10 and 11 (W = 1.1 nm and 1.3 nm) for armchair 
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nanoflakes and N—6 (W—1.3 nm) for zigzag nanoflakes. 
To keep the computational burden low, we chose few rep- 
resentative lengths (L ~ 3.2 nm to 6.1 nm). 



for the fact that the bilayer orbitals play a role in the 
localization/delocalization process and the degree of lo- 
calization will get affected due to interlayer coupling. 
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FIG. 5: (Color online) Variation of the armchair and zigzag 
nanoflake gaps with different lengths are shown. In longer 
nanoflakes, the non-magnetic gaps are shown to be smaller 
than their magnetic counterpart. 

It is seen from Fig. [3] that non- metallic gaps are weakly 
dependent on the nanoflake lengths whereas the metallic 
nanoflake gaps show strong length dependence. More- 
over, it is clearly seen that interlayer magnetism enhances 
the intrinsic gaps in longer nanoflakes. In Ref. 10, it 
was suggested that in non-metallic monolayer armchair 
nanoflakes, the 7r-orbitals are strictly localized along the 
zigzag edges and therefore the gaps are insensitive to the 
increase in length, whereas for metallic nanoflakes, the 
orbitals are delocalized throughout the flake and there- 
fore are sensitive to changes in length. We believe that 
this will also be true for the bilayer nanoflakes except 



C. Electric field effects on the gaps 

In this section, we discuss the effect of external elec- 
tric fields, applied perpendicular to the layers, on the 
gaps of armchair and zigzag nanoflakes. We considered 
four different values of the electric fields below the SiC>2 
dielectric breakdown field of 1 V/nm. Both wide and nar- 
row nanoflakes were chosen with L « 3.2 nm for armchair 
and L « 3.3 nm for zigzag nanoflakes. The electric field 
decreases the gap in both armchair and zigzag nanoflakes 
with a-alignments (Fig. a behavior we predicted for 
bilayer nanoribbons.— A similar behavior is seen for /3- 
aligned nanoflakes (figures not shown). 
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FIG. 6: (Color online) Variation of the energy gap with 
perpendicular external electric field for bilayer nanoflakes 
with a-alignments for (a) armchair (L ~ 3.2 nm) and (b) 
zigzag nanoflakes (L ~ 3.3 nm). For each class of armchair 
nanoflakes, three representative widths specified by the in- 
terger p=3, 6, and 8, are chosen. The gaps are shown at three 
representative values of the external electric field strength. 
The maximum electric field strength applied is close to 1 
V/nm. 

We note that, recently, we predicted the existence of a 
critical gap of about 0.2 eV for bilayer nanoribbons below 
(above) which the electric field has the effect of increasing 
(decreasing) the gap^ We could not verify that such a 
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critical gap can exist for armchair and zigzag nanofiakcs 
because of gap values, all exceeding 0.2 eV, for the chosen 
widths and lengths. We did additional calculations of 
electric field effects on the longer nanoflakes up to £=6.1 
nm to search whether there exists a critical gap but again 
due to the large gap values in longer nanoflakes, no such 
critical gap was found. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied, using a first principles 
DFT method, the tunability of bilayer nanoflake gaps as 
a function of interlayer magnetism, lengths and widths, 
and external electric fields. Bernal (or AB) type of inter- 
layer stacking with two edge types (armchair and zigzag) 
and two edge alignments (a and (3) were considered. We 
identify three distinct classes of armchair gaps and show 
that interlayer magnetism plays an important role in en- 
hancing the confinement-induced gaps. Length-confined 
metallic armchair nanoribbon gaps are strongly affected 
by variations in the length. The energy gap as a function 
of the applied electric field show decreasing trend for both 



the edge types and alignments. However, the existence 
of a critical gap can not be ruled out for nanofiakcs with 
very small intrinsic gaps. We expect that the present re- 
sults will help stimulate further studies of bilayer gaps 
in the presence of additional external perturbations such 
as a substrate and contacts, and motivate experiments 
to unravel the complicated interplay of magnetism, geo- 
metrical confinement and edge type. 
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